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Abstract. Electromagnetic superradiant field coherence exists in a condensed 
matter system if the electromagnetic field oscillators undergo a mean displacement. 
Transitions into thermal states with ordered superradiant phases have been shown 
to theoretically exist in Dicke-Preparata models. The theoretical validity of these 
models for condensed matter has been called into question due to non-relativistic 
diamagnetic terms in the electronic Hamiltonian. The microscopic bases of Dicke- 
Preparata thermal superradiance for realistic macroscopic systems are explored in this 
work. The impossibility of diaelectric correlations in condensed matter systems (via 
the Landau-Lifshitz theorem) provides a strong theoretical basis for understanding the 
physical reality of condensed matter thermodynamic superradiant phases. 



PACS numbers: 78.60.Kn, 78.60.Kn.Fi, 78.70.-g 



1. Introduction 

In quantum electrodynamic (QED) theory, the magnetic field B and the Maxwell 
displacement field D are non-commuting quantum fields; i.e. in Gaussian units 

[A(r), B,{r')] = -Amhce,,A6{v - r'). (1) 

The electromagnetic field in a condensed matter system is said to be coherent if and 
only if (on average) B 7^ or D 7^ or both situations hold true. In terms of the 
electromagnetic oscillator modes, there is QED coherence when one or more of the 
long wavelength field oscillators exhibit a spontaneous non-zero mean displacement (even 
in the absence externally applied fields). Ferromagnets and ferroelectrics clearly exhibit 
such spontaneous electromagnetic field ordering on a macroscopic scale. Microscopic 
QED coherence may also occur on the smaller length scales of atomic and molecular 
physics as well as on the length scale of mesoscopic coherent domains. For example, 
atomic or molecular magnetic moments clearly exhibit finite magnetic fields B(r) 7^ 
in a spatial domain surrounding the moments. Dicke models of thermal superradiance 
have been extensively studied [2-23] and constitute an important class of ordered 
electromagnetic field coherent states. 
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The Dicke Hamiltonian [p^ for a single photon oscillator may be written in the form 

Hoicke = 2 (-^^ + ^ooQ'^) ~ fQ + Helectronic- (2) 

The first term on the right hand side of Eq.(|^) represents a single photon oscillator 
mode. The electronic states are described by the Hamiltonian Heiectronic- In the Dicke 
model, Heiectronic describes a set of two (energy) level molecules. The oscillator force 
operator / is determined by the electronic dipole moment operators of the molecules. 
The existence of a superradiant phase transition has been rigorously proven ^ 
for the strong coupling Dicke Hamiltonian in Eq.(H). The reason for the superradiant 
phase transition may be understood by computing the shift in the photon oscillator 
frequency from ojoo to ojq due to the dipole interaction. The renormalized frequency ojq 



is determined by the oscillator strength sum rule ||25 



2 

Q + - T{uj)duj where T{uj) > 0. (3) 
Jo 

The oscillator frequency dependent damping function, denoted by T^u) and caused by 
the random force /, lowers the square of the frequency from ut'^ to Uq. If the damping 
is sufficiently strong, then uJq < and the oscillator becomes unstable about the old 
equilibrium position Qoid = 0. A new equilibrium position is reached Qnew 7^ which 
describes the superradiant ordered state. 

There is little doubt that strong damping can produce a superradiant ordered 
phase if the Dicke Hamiltonian is employed. However, the validity of the Dicke model 
Hamiltonian as a description of real condensed matter systems has been questioned [27- 
31]. At issue are the diamagnetic terms in a more realistic microscopic non-relativistic 
electronic Hamiltonian [32-35]. With B = curl A, there are in the non-relativistic 
electronic Hamiltonian both linear and quadratic terms in B. The quadratic terms 
in B induce possible diamagnetic effects in condensed matter systems. Diamagnetism 
is dominant in some ordered phases of matter (e.g. superconductors) but ferromagnetic 
correlations are dominant in some other ordered phases of matter (e.g. ferrites). For 
magnetic field coherence B 7^ 0, the frequency shift in Eq.(^ is renormalized (in part) by 
diamagnetic strengths independent of the oscillator damping coefficient T{uj). However, 
the superradiant ordered phases in the Maxwell displacement field D are not adversely 
affected by diamagnetic couplings. 

The essential physical ideas have been explored by Landau and Lifshitz The 
conventional QED Hamiltonian contains only linear interaction terms in D. Landau and 
Lifshitz prove the following important theorem for materials described by a dielectric 
constant e (wherein D = eE) and a magnetic permeability fi (wherein B = /iH): 

Theorem 1 The magnetic permeability /i > can be either paramagnetic fi > 1 or 
diamagnetic n < 1. The dielectric constant must be paraelectric e > 1. There is 
no "diaelectricity" in condensed matter systems, i.e. the range 1 > e > Q is strictly 
forbidden. 
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The above theorem and the generahzations to be discussed below provide a strong 
theoretical basis for physical reality of the thermal superradiant phase. Our purpose 
is to discuss the theoretical microscopic bases of Dicke-Preparata models as applied to 
realistic macroscopic condensed matter systems. 

In Sec. 2 we describe the conventional QED Hamiltonian. The electronic degrees of 
freedom are considered to be non-relativistic. A canonical transformation is introduced 
which allows for a precise microscopic definition for the transverse Maxwell displacement 
field D = E + AttP. In Sec. 3 the notion of diabatic and adiabatic changes in 
electro dynamic processes will be defined. Dissipation is introduced as a diabatic (i.e. 
non-adiabatic) process. In Sec. 4, the statistical thermodynamics of the electric dipole 
moment interactions will be discussed and a generalization of the Landau-Lifshitz 
Theorem |l| will be proved. The diabatic dissipation will be related to the thermodynamic 
dielectric response of the medium. In Sec. 5, it will be shown how the instability of 
the thermodynamic dielectric response function is the signature of a transition into a 
superradiant phase. 



Hqed 



2. QED Hamiltonian 

If we choose a vector potential in the Coulomb gauge, 

B = curlA, divA = 0, (4) 
then the QED Hamiltonian of interest in the work which follows has the form 

1 1 {\Ef + \B\^) dh + H[A], (5) 

where E' is the operator which denotes the transverse part of the electric field, 

E» = 4,,fc (^) = 4.>fc curl (^) , (6) 

and H[A] denotes the Hamiltonian (including Coulomb interactions) of the charged 
particles. The Coulomb Hamiltonian when B = has the form 

ff|A = 01 = i/cw = -E|^A,-5:^A„ + C/, (7) 

j a. 

where the total potential energy for the electrons and nuclei in the condensed matter 
system is given by 



Kj<k •' a<b ja ■' 



U = e^ 

For A ^ 0, 
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where Aj = A(rj), = A(Rq) and the interaction of the magnetic field with the 
particle spins is given by 



\2mc/ 



^s,-B(r,)-5^ 



gaZa\e\ 

2M„c 



Sa ■ B(Ra). 



(10) 



Let us consider a basis in which the Coulomb Hamiltonian is diagonal; i.e. 

H[A = 0]ij, = Hcouii^, = Wj[A = 0]^j where j = 0, 1, 2, . . . . (11) 
In the Coulomb basis, one may define the matrix elements 

H,k[A] = {^k,H[A]^,) (12) 
so that the Hamiltonian in Eq.(^ has the equivalent matrix representation 

fHoo[A] iJoi[A] Ho2[A] ...\ 
H^o[A] Hu[A] H^2[A] ... 
H2o[A] H2i[A] H22[A] .., 

V : ! I ••./ 

In principle, the Hamiltonian can be brought to diagonal form by a unitary 
t r ansf ormation 

W[A] = U^A]H[A]U[A] where U^[A] = [/"^[A] (14) 



H[A] 



(13) 



and 



W\A] 



fWo[A] 
Wi[A\ 
W2[A] ... 



(15) 



In virtue of gauge invariance 

W[A] = W[B]. (16) 

The unitary transformation can be employed to transform the total Hamiltonian in 
Eq.(^ into the adiabatic representation 

n = U^[A]HQEDU[A]. (17) 

When acting on the electric field in Eq.(^, the transformation defines the Maxwell 
displacement field D and polarization P via 

U^[A]E'U[A] = E = D - 47rP, (18) 

where 



Dfr) 



A-Kihc 



6 



and 



P(r) = -ihcU^lA 



5A{r 

fSU[A]\ 



Anihc curl 



\6A{r) J 



ihc 



5B(r)7' 
SW[A 



6A{r) 



U[A]. 



(19) 



(20) 
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Eqs.(||), and (JT^) - ( pOD imply the adiabatic Hamiltonian representation 

'H = ^J (|D-47rPp + |Bp)ci^r + iy[A], (21) 
which may be conveniently written as 

n = ^ J |D|2rf3r - y D • Pdh + H'[A]. (22) 

In Eq.(|2g), 

H'[A] = ^ J \B\'^dh + W[A] + 2TT J \P\'^dh. (23) 

Finally, the transverse current operator in the adiabatic representation of Eqs.(14) and 
(p!5|) is given by 

J(r) = -cUH^]^-^U[^]- (24) 

The notion of "diabatic damping" is associated with a decomposition of the current 
Eq.(p^ into adiabatic and non-adiabatic parts. 



3. Damping and Electrical Conductivity 

For the unitary transformation in Eq.(|l^) we have the differential identity 

U^^-^U = 4r iU^HU) - U^HUU^^-^ - ^-^UU^HU, (25) 
which together with Eqs.(0) and (|2^) yields the current decomposition 

= ^^^^ 

One may employ the gauge invariant B = curl A in the form [see Eq.(|16])] of an 
adiabatic current 

5W[A] ,(5W[B] 
•J adiabatiA'^) = — c ^ . , , = — c curl——— = c curlM. 27 

The gauge invariant adiabatic magnetization has a form which follows from Eqs.(|T^) 
and (13); i.e. 

f6Wo[B]/6B 

5iyi[B]/5B 

6W2[B]/6B ... 

V : ; : ••/ 

The dissipative (or "diabatic") part of the current operator is given by 

Jd(r) = ^[l^,P(r)] = P(r). (29) 

Eq.(p6D then takes the conventional form having both magnetization and polarization 
parts; i.e. 

J(r) = c curlM{r) + P(r). (30) 



M = - 



(2^ 
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Although Eq.(|30|) is well known in classical electrodynamics, we have given the proof 
from a fully quantum electrodynamic viewpoint. 

The diabatic part of the current in Eq.(^) describes the dissipation via the 
non-local transverse electrical conductivity tensor (Tjj(r, r', + iO"*"). The microscopic 
expression for the transverse conductivity is determined by the fluctuation dissipation 
theorem 

1 j-P POO 

a,,(r,r',C) = T / / e*^* ( Jrf,,(r', -zA) Jrf,,(r, t)) rftdA. (31) 
" Jo Jo 

Note: (i) The complex frequency obeys 53m C > 0. (ii) The thermal average (...) is over 
the charged particle degrees of freedom, (iii) The time variation of operators in Eq. ( PT| ) 
employs the Hamiltonian -ff'[A] in Eq. (P^D . (iv) Finally, 



For the transverse dielectric properties of charged particles, one employs the electric 
susceptibility 

• /.oo 

X^J{r,r'X) = -J^ e<*([P,(r,t),P,(r',0)])dt. (33) 
Integrating Eq.(|33|) by parts yields 

- ^Cx^J{r, r', C) = h,,{r, r') + ^ e^^* ( [^(r, t), P,(r', 0)] ) dt, (34) 
where the equal time commutator contribution is given by 

h.Ary)='-mr),P,{r')]). (35) 
Employing the Kubo-Martin-Schwinger [37-40] condition 

t ( [P,(r, t), P,(r', 0)] ) = J^^ (P,(r', -?A)P,(r, t)) d\ (36) 

and Eq.(|25|) in Eq.(^) yields a simple relationship between the conductivity and the 
susceptibility 

- <Xii(r, r', C) = aij{r, r', C) + h,j{r, r'). (37) 
The dissipative part of the conductivity is determined by 

^e{cr,j{r,r',uj + iO+)} = ujQm{xij{r,r' ,uj + iO+)}. (38) 

The quantum mechanical hall conductivity contribution hij{r,r') from the transverse 
polarization is determined by the commutation relation in Eq. (p5|) . 

Note that the response functions discussed above are defined with respect to an 
applied magnetic field B = curlA. Magneto-conductivity in aij and the Faraday effect 
in Xij are implicitly included in the above considerations. For superradiance described 
by a coherent field D 7^ one averages over magnetic field fluctuations. The mean 
magnetic field (on macroscopic length scales) obeys B = 0. Let us now consider the 
details of the magnetic field fluctuations. 
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4. Statistical Thermodynamics 

The quantum electrodynamic free energy F of a condensed matter system can be written 
as 

g = -kBTln {Tr e-^/^'^^j (39) 

where H is given in Eq. (p2|) . The complete trace over quantum states includes 
electromagnetic field degrees of freedom as well as charged particle degrees of freedom. 
Thus 

Tr e-^^l^-^ = Trp,A) {Tr^cMr.e,) e'^/'^-^) . (40) 

If one treats the electromagnetic field trace in the quasi-classical limit of a functional 
integral 

/"tt fVDMV A(r)\ , , 

r ^ / 

then Eqs.(P^ and (^)-(^) imply the functional representation 

-"(-^j-Z-P^^)!!™. (42, 

The total free energy functional of the Maxwell displacement field obeys 

^to^D] = ^ [ \Bfd\ + ^[D], (43) 



Stt 



where 



and 



H"\A, D] = H'\K\ - j P(r, [A]) ■ B{r)d^r. (45) 



The coupling into the Maxwell displacement field in the Hamiltonian Eq. (|45|) is 
linear. Eqs.(^^ and (|45|) imply that the free energy J^lD] is a convex upward functional. 
For the second functional derivative 

PI) 

we have the following: 

Theorem 2 The spectrum of the zero frequency susceptibility {xa[D]} defined by the 
eigenvalue equation 

J2 I X.,(r,r', [D])ei(r')ciV = XA[D]a(r) 
obey Xa[D] > for all A. 
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In order to prove a generalization of the Landau and Lifshitz Theorem |1|, one must 
investigate the second functional derivatives 



= A,,(r-r')-47rx.,(r,r',[D]), (47) 
where Eq.(^) has been invoked and the transverse delta function is defined as 



3 



The spectrum of the zero frequency transverse dielectric constants {£:a[D]} is defined 
by the eigenvalue equation 



J2 J ^..(r,r',[D])ei(r')rfV= (-i^) ^IW 



and obeys 



(49) 



(50) 



_l-47rxA[D] 

The thermodynamic stability condition for the free energy jFtot[D] can be 
summarized by the following: 

Theorem 3 Thermodynamic stability for a Maxwell displacement field D requires that 
ex[B] > for all A. 

Proof: For a Maxwell field D to represent thermal equilibrium, one expects the free 
energy jF^ifD] to be at a minimum. The second (functional) derivative conditions for 
achieving a thermodynamic free energy minimum are those stated in the theorem in 
virtue of Eqs.(0) and (||). 

The central result of this section is the generalization of the Landau-Lifshitz 
Theorem 

Theorem 4 Thermodynamic stability for a Maxwell displacement field D requires the 
paraelectric inequality e\[D] > 1 for all A. Diaelectric behavior (0 < e\\D] < 1) is 
strictly forbidden. 

Proof: Theorem |^ follows directly from Theorem Eq.(|50D and Theorem 0. 

We now reconsider the dissipative properties of a condensed matter system after 
coherent averaging over magnetic field B = curlA fluctuations. We employ the 
functional averaging measure /(...) Hrl^^/^^^'^) Eqs.(P3D and (p^. For 

example, the new susceptibility 

X.,(r,r',C,[D]) = ly^ e<*(([P.(r,t),P,(r',0)]))rft, (51) 

where the complete "double averaging" ((...)) is over both charged particle motions and 
magnetic field fluctuations with the Maxwell displacement field D held fixed. The time 
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dependent operators in Eq.(|5lD are with respect to the Hamiltonian H"[A,Y)] defined 
in Eg. (|i5|) . The dissipative motions of the condensed matter system are then described 
by 

^e{a,j{r, r', cu + zO+, [D])} = cj^>m{xii(r, r', to + iO+, [D])}. (52) 
Finally, the static susceptibility in Eq.(^B|) is the zero frequency limit 

X.,(r,r',[D]) = lim3?e{x.,(r,r',cu + zO+,[D])}. (53) 

The thermodynamic stability test for a coherent electromagnetic superradiant state 
follows from Theorem Eq. (|50|) and Theorem ^ 



5. Thermodynamic Stability 

In the absence of external electric fields, the stable values of D are such as to minimize 



36| the total free energy ^^^[D]. The free energy minimization condition implies a 



vanishing thermal electric field; 



The first derivative Eq.(|5^) will in general have more than one possible solution with a 
non-zero Maxwell displacement field [D] . The multiple solutions correspond to differing 
possibilities for superradiant coherent domains. A necessary condition for a true free 
energy minimum has been proved in Theorem ^. Employing Eq. (|50|) we find the stability 
condition 

< 47rxA[D] < 1. (55) 
From (i)Eq.(^ll), (ii) the resulting dispersion relation 

Xij(r,r ,C, [D]) - ' 



Xij{r, r', [D]) = lim Xii(r, r', C, [D]) 

2 1"°° duj 
= -/ 53mx.,(r,r',a; + zO+,[D]) — , (56) 

(iii) Theorem^ and (iv) the stability Eq.(^), it follows that strong dissipation, i.e. 
substantial {^mx{oj + iO^)/^^}, tends to yield thermodynamic instabilities. 

Symmetry under parity transformations yield the solution D = corresponding to a 
"normal" phase. If 47rx;v[D = 0] > 1 for some A, then the normal phase is unstable. The 
true free energy minimum will arise for coherent superradiant domains with [D ^ 0]. 
There may be many possible superradiant domain configurations as discussed above. 
Domain walls and/or normal phase regions may separate the superradiant domains 
wherein D = 47r ((P(r))) ^ 0. 
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We have discussed a generalization of the theorem by Landau and Lifshitz whereby e > 1 
for materials described by D = eE. The theorem and its generalization is independent 
of magnetic permeability /x when the material also obeys B = /iH. Magnetic properties 
may be either paramagnetic or diamagnetic, but diaelectric properties of matter are 
ruled out by quantum statistical mechanical stability considerations. 

The stability criteria are crucial for the understanding of stable superradiant 
domains. In particular, for all of the possible eigen modes of the dielectric susceptibility, 
the generalized Landau-Lifshitz stability condition £a[D] > 1 is crucial for testing 
whether (or not) a superradiant domain is thermodynamically stable. Contrary to what 
has some time ago appeared in the literature p9| the dielectric stabihty conditions are 
independent of possible diamagnetism. No approximations to the free energy functional 
JFjo4[D] have been required to reach important general conclusions, however further 
approximations are required for concrete computation. 
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